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Phase-separation fronts leave in their wakes morphologies that are substantially different from 
the morphologies formed in homogeneous phase-separation. In this paper we focus on fronts in 
binary mixtures that are enslaved phase-separation fronts, i.e. fronts that follow in the wake of a 
control-parameter front. In the one-dimensional case, which is the focus of this paper, the formed 
morphology is deceptively simple: alternating domains of a regular size. However, determining the 
size of these domains as a function of the front speed and other system parameters is a non-trivial 
problem. We present an analytical solution for the case where no material is deposited ahead of 
the front and numerical solutions and scaling arguments for more general cases. Through these 
enslaved phase-separation fronts large domains can be formed that are practically unattainable in 
homogeneous one-dimensional phase-separation. 



I. INTRODUCTION 

Phase-separation fronts are critically important for 
the formation of phase-separation morphologies in many 
practical applications. Two industrial examples are the 
domain formation in binary alloys [l[ and immersion pre- 
cipitation membranes @, M 0| • A binary iron-carbon 
alloy which is slowly cooled from one side of the sam- 
ple below its critical temperature may form a structure 
of alternating carbon-rich and carbon-poor bands 0, Q . 
The temperature change moving through the alloy in- 
duces phase-separation where it crosses a critical value. 
A very different example is the formation of immersion- 
precipitation membranes. Here a polymer-solvent mix- 
ture is immersed in a non-solvent. Phase-separation is 
induced as the solvent diffuses out of the mixture and is 
replaced by the non-solvent. Phase-separation starts at 
the initial interface. From there a phase-separation front 
advances into the mixture. Both of these examples have 
in common that a space and time dependant parame- 
ter, i.e. the temperature in the binary alloy case and 
the solvent concentration in the immersion precipitation 
case, induce phase-separation. We call such a quantity 
a control parameter. The advance of the temperature 
or the solvent concentration then constitutes a control 
parameter front which is responsible for inducing phase- 
separation. 

Despite its critical roll in the formation of complex 
structures, only limited attention has been paid to this 
phenomenon. There are a few notable exceptions. Fu- 
rukawa [7[ examined the phase-separation in two dimen- 
sional binary mixtures following a sharp control param- 
eter front and noted that the morphology formed was 
strongly determined by the speed of the front, and that 
morphology may change quite suddenly if the front speed 
is changed. Hantz Q rediscovered the front speed de- 
pendence of morphology using a novel rotating control 
parameter front where the front speed depends on dis- 
tance from the axis of rotation. Dziarmaga suggested 
a model that is very similar to the one presented in the 
current work as a model for phase-separation in the early 
universe, and found that for slow fronts the density of do- 



main walls is linear with front speed, i.e. domain size is 
proportional to the inverse of front speed. 

When a control parameter front is not progressing too 
quickly a phase-separation front follows directly behind 
and therefore moves at the same speed. We call those 
fronts enslaved phase-separation fronts. In the opposite 
case, the phase-separation front will not keep up with 
the control parameter front. This then results in a free 
phase-separation front moving into an unstable medium. 
This special case has received a significant amount of 
attention and a good review is given by Saarloos [i"ol |. 

Enslaved phase-separation fronts give rise to a host 
of complicated phenomena, as witnessed by the large 
variety of morphologies formed. For example, those 
formed by the previously mentioned immersion precipita- 
tion membranes [2| . These morphologies are found to de- 
pend strongly on processing conditions which correspond 
to different shapes and speeds of the control parameter 
fronts. This is why we initially focus on the simplest 
example of enslaved phase-separation fronts. Here we 
consider only one-dimensional abrupt (i.e. not spatially 
extended) control parameter fronts in binary mixtures. 
This already gives rise to unexpected and non-trivial be- 
havior as explained below. In particular it is noteworthy 
that large domains can be formed efficiently, even though 
such domain sizes are practically unobtainable through 
homogeneous phase-separation in one-dimensional sys- 
tems. 

In the next section we review the theory of spinodal de- 
composition, as it is relevant to the understanding of the 
dynamics of phase-separation at the front. We then de- 
fine our model which consists of a Cahn-Hilliard equation 
with an underlying 4 free energy. The control parame- 
ter front enters though the time and space dependence of 
the parameters in the free energy. We expect the details 
of the phase-separation dynamics to depend strongly on 
the shape of the control parameter profile and its dynam- 
ics. To simplify our analysis in this paper we focus on 
the simplest case consisting of a control parameter front 
which is a sharp step moving with a constant velocity. 
To numerically explore the dynamics of this model we 
developed a lattice Boltzmann method (LBM) simula- 
tion. Generic results of these simulations elucidate the 
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general behavior of an enslaved phase-separation front. 
As the control parameter front advances into the mixed 
region, phase-separation ensues. Excess material is then 
diffused both ahead of the front and through the form- 
ing domain into the next domain. This continues until 
switching occurs and a new domain of the opposite kind 
is formed. This is a regular process that repeats, result- 
ing in regular domains after a few cycles. 

To predict the morphology we only need to find the 
dependence of this one length-scale as a function of the 
nondimcnsional parameters of the equation of motion. 
Despite its apparent simplicity this is still a non-trivial 
problem. We show that this can be done analytically in 
the case where the mobility ahead of the front is negligi- 
ble and the front is moving very slowly. The main ingre- 
dients in the analytical solution are solving for the time 
dependence of the concentration right behind the front 
and the condition for the switching to a new domain. 
Together, those allow us to determine the switching time 
and hence the domain size. 

Using the results from this special case as a founda- 
tion, we examined more general enslaved fronts. When 
we investigated the effect of non-zero mobility ahead of 
the front we observed that the domain-size has an unex- 
pected non-monotonic dependence on the mobility. We 
were able to qualitatively explain this behavior as a com- 
petition of the diffusive deposition of material ahead of 
the front and the switching condition. 

With satisfactory results for this model in one dimen- 
sion, we conclude with a discussion of possible extensions 
such as using the front speed to imprint a domain pat- 
tern, and a roadmap of future research of enslaved front 
systems. 



II. PHASE-SEPARATION AND FRONTS 

As previously mentioned, phase-separation occurs af- 
ter some physical control parameter, such as tempera- 
ture, changes so that the system becomes thermodynam- 
ically unstable. The system forms two or more distinct 
coexisting phases, and the property which distinguishes 
the phases is called the order parameter. We are con- 
cerned with systems that transition from a one-phase 
state to a two-phase state as a control parameter is var- 
ied. Despite the fact that most theoretical work on phase- 
separation has focused on the homogeneous case, most 
practical cases exhibit phase-separation that does not oc- 
cur everywhere-at-once, but starts at one or more initial 
points and spreads through the system. The boundary 
between phase separated and non phase separated mate- 
rial is called the phase-separation front. Likewise, when 
the control parameter changes inhomogeneously we refer 
to the boundary between the single phase and two-phase 
regions as the control parameter front. The correlation 
of the phase-separation front and the control parameter 
front depends on the exact nature of the dynamics of the 
control parameter. 



Consider a binary mixture of ,4-type and B-type ma- 
terial described by the free energy: 



F 



dx 



(1) 



This is a general fourth-order expansion of more specific 
free energy expressions. The cubic term which would 
normally appear has been scaled away, without loss of 
generality, by choosing the critical concentration to be 
zero. The c parameter, which can usually be neglected, 
can matter here because it may take different values on 
either side of the control parameter front. The dynamics 
of phase-separation is diffusive and is described by the 
Cahn-Hilliard equation 



d t 4> = V(mVfi) , 



(2) 



where the chemical potential is derived from the free en- 
ergy as 



dF 

p = — = a<j) + b<j) 3 - kV 2 + c . 



(3) 



Here a is the control parameter that determines the sta- 
bility of the system; a > corresponds to the one-phase 
region and a < corresponds to the two-phase region. 
The order parameter <f> = pj± — Pb is the concentration of 
the two component materials. For homogeneous param- 
eters the linear c term can be neglected because it does 
not change the equilibrium properties or the dynamics 
of homogeneous systems. However, when c = c{x) it will 
introduce chemical potential gradients between the mate- 
rial on either side of the control parameter front. Of the 
other parameters, m is the mobility, k determines the in- 
terfacial energy cost, and the amplitude of the non- linear 
term b determines the equilibrium values of the order pa- 
rameter: 



(4) 



We now discuss different phase-separation mechanisms 
described by this model. 



A. Homogeneous Control Parameter Change 

If the control parameter does change rapidly every- 
where at once — the well studied homogeneous quench — 
no control parameter front exists. If the system is also ho- 
mogeneous in composition a phase-separation front will 
not form, in which case phase-separation occurs via spin- 
odal decomposition or nucleation. 

Another interesting situation occurs when the control 
parameter is changed homogeneously but gradually. In 
such a situation, cascades of sequential phase-separation 
events are observed, as shown in the work of Vollmer et 
al. [ll], Q3 • Such phenomena are probably also exhibited 



3 



for extremely extended enslaved control parameter fronts 
but are outside the scope of the current paper. 

In the case of a homogeneous quench, typically spin- 
odal decomposition is observed. It manifests itself in an 
exponential change of the order-parameter from initial 
fluctuations towards one of the two equilibrium values. 
The growth rate depends on the fluctuation's wavelength, 
and the fastest growing wavelength will outgrow the oth- 
ers and quickly dominate the morphology. This fastest 
growing wavelength is called the spinodal wavelength, de- 
noted A sp . Though it is not often discussed in the anal- 
ysis of homogeneous quenches, we refer to the reciprocal 
of the fastest growth rate as the spinodal time, denoted 

To derive the spinodal wavelength and time one per- 
forms a linear expansion of the Cahn-Hilliard equation 
(|2j) around the initial concentration <pi n , then Fourier 
transforms this linearized Cahn-Hilliard equation into k- 
space which results in the ODE: 

d t 4>{k) = -m [(a + Zb<p 2 n ) k 2 + nk A ] 4>(k) = R(k)4>(k) . 

(5) 

The solution of this equation is the exponential growth 
of fluctuations where the growth rate depends on the 
angular wavenumber k: 



ftk) 



(6) 



The angular wavenumber with the fastest growth rate 
corresponds to the spinodal wavelength: 



A., 




(7) 



The fastest growth rate corresponds to the spinodal time 
1 4k 



R(k sp ) m(a + 3&</>; 



2 ^2 



(8) 



Such a morphology, with an initial domain wavelength of 
A sp , will coarsen to larger length-scales, but the process 
is inefficient in one dimension since it obeys a logarithmic 
growth law L oc \n(t) at late times Therefore it is 

hard to form large structures in a finite time. 

From the spinodal wavelength and time we define a 
spinodal speed 



*s P irm , , 
= — - = —= [—a — 3b(p 



.3/2 



2 k 



(9) 



which can be thought of as a natural speed of phase- 
separation. This is important to our later analysis of the 
dynamics of enslaved phase-separation fronts because it 
allows us to define a nondimensional front speed. 



much more rapidly than through spinodal decomposition. 
From this defect a front spreads through the system. Sar- 
loos et al. performed an analysis of the speed of such 
fronts for many different types of systems. For a con- 
served order parameter system described by the Cahn- 
Hilliard equation J5]| with a — —I and m = b = k = 1, 
the speed of the front, which they call the linear spread- 
ing speed, is: 



' 34 + UV7 
27 



(1 - 3#„) 



3/2 



0.73 u s 



(10) 



The initial domain wavelength is found to be smaller than 
the spinodal wavelength A* « 0.35 A sp and the morphol- 
ogy then coarsens in time [Iol |. 



C. Enslaved Phase-Separation Fronts 

Freely propagating phase-separation fronts can also 
form as a special case in a system with a control parame- 
ter front. If the control parameter front moves at a speed 
u > u*, this can result in the suppression of spontaneous 
spinodal decomposition, yet not inhibit the advancement 
of a freely propagating phase-separation front. However, 
if the control-parameter is moving slower than the linear 
spreading speed u < u* , the phase-separation front can- 
not propagate freely and becomes enslaved by the control 
parameter front. 

The nature of a control parameter front highly depends 
on the particulars of the physical system. For instance, 
if temperature is the control parameter of a long binary 
metal alloy rod which is being cooled from one end, the 
control parameter front will resemble the familiar error- 
function solution of the heat-diffusion equation. How- 
ever, if the same rod is being extruded from a hot oven 
into a cooling environment, the control parameter front 
will be a much sharper transition. The shape and speed 
of the control parameter front will have an impact on the 
morphology formed in the wake of the front. In this pa- 
per we focus on phase-separation fronts enslaved by an 
abrupt control parameter front moving with a constant 
velocity. This abrupt control parameter front will be rep- 
resented by a step function for the parameters in the free 
energy and the mobility in the Cahn Hilliard equation. 

We now introduce a lattice Boltzmann method to 
numerically investigate this system. We will see that 
this control-parameter front induces phase-separation dy- 
namics which is confined to a limited region around the 
front. Consequently we can reduce the computational 
cost by focusing on a narrow region around the front. 



B. Free Propagation of a Phase-Separation Front 

Phase-separation fronts can form if a system has be- 
come unstable due to a sudden homogeneous control pa- 
rameter change and a front is nucleated. This can oc- 
cur if a defect causes local phase-separation to occur 



III. LATTICE BOLTZMANN METHOD FOR 
PHASE-SEPARATION FRONTS IN ID 

As we mentioned above, the dynamically important re- 
'ion for this model is restricted to a narrow area around 
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the front. So rather than simulate a stationary system 
with a front moving through it, we consider the equiv- 
alent case of material advected with velocity u past a 
stationary phase-separation front. This allows us to use 
a much smaller simulation size. This changes the dif- 
fusive equation of motion ([2|) to the advection-diffusion 
equation 



dt<j> + V(ucf>) = V(mVju) , 



(11) 



where u was the constant phase-separation front speed, 
which now appears as the material advection speed. The 
parameters which define the abrupt control parameter 
front are 

a(x) = a s + (a M - a s )Q(x - x ) , (12) 
b(x) = b s + (b M - b s )Q{x - xo) , (13) 
m(x) — m s + (m,M — ms)0{x — x ) , (14) 
k(x) = k s + (k m - ks)Q(x - x ) , (15) 

where Q(x) is the Heaviside step function and xq is the 
position of the front. Our naming convention uses the 
subscript M to denote the mixing region ahead of the 
control parameter front, and the subscript S to denote 
the separating region behind the control parameter front. 

The LBM implementation of the advection-diffusion 
equation is similar to that of the ordinary diffusion equa- 
tion, however the key differences of an advecting mate- 
rial with spatially dependent parameters warrents further 
clarification. 

The standard lattice Boltzmann equation with the 
BGK approximation is 

fi(x + v i ,t + l)-fi(x,t) = ±[f?-f i (x,t)] , (16) 

for discrete integer time, and a finite set of discrete ve- 
locities Vi defined such that x and x + Vi are sites on 
the spatial lattice. To determine the macroscopic evo- 
lution equation for this model we define a continuous 
fi = fi(x,t), and expand the first term to second order 
in a Taylor series as 



fi(x+v t ,t + i) 

where D is the total derivative operator: 



fi + Dfi + -D 2 f 



0(D 3 



D = d t + v ia V a . 
This allows us to rewrite Eq. (|16p as 

Df i + ±D*f i + 0(D s ) = ±{f° 
and to first order we obtain: 

fi = /? - 7 I) I, 



0(D 2 



(17) 



(18) 



(19) 



With repeated use of Eq. to replace f in Eq. lfT8]) 
with the local equilibrium functions. The hydrodynamic 
limit of the lattice Boltzmann equation then becomes 



D/r- D[r-\)Df): 



0(D 3 ) = -(/?- /;) , (20) 



with the implicit generalization that that r = t(x, t). We 
then define the equilibrium moments 

I> = E/° = 0, (21) 

i i 

J^fSa = SU a (f>, (22) 

i 

y j fiViaVif] = Sfi + S 2 U a U0(f) , (23) 



corresponding to a conserved order parameter and a cur- 
rent of su(j>. Summing up the Taylor expanded LBE (|20p 
over i we obtain the equation of motion 



_1 3t0 + V(u0) + 0(d 3 ) = V(mV/i) 



(24) 



This is Eq. pd|) to second order, with the addition of a pa- 
rameter s which allows for rescaling simulation time. The 
ability to scale simulation time allows us to trade compu- 
tational speed for enhanced numerical stability and im- 
proved Galilean invariance [l4j. In this formulation the 
mobility m is given by: 



m(x) = t(x) - - . 
In one dimension we use three discreet velocities 

t>i = {0,-l,+l}, 
which define the set of equilibrium distributions 
/ ° = (1-sV^-s/i, 
/-i = \{{s 2 u 2 ~su)^ + S ^) , 

/+i = i((sV + «0^ + «/i) • 



These equilibrium moments are then used with Eq. (fTf 
to calculate the time evolution of the system. 

We now have to pay the cost for smaller simulation size 
allowed by fixing the position of the front in terms of more 
complicated inflow and outflow boundary conditions. To 
represent a front moving in the positive x direction we 
define u < which moves the material from right to 
left past the stationary front. This defines a fixed inflow 
boundary at x = x max and free outflow boundary at 
x = 0. Our inflow boundary condition is then given by 



(25) 

(26) 

(27) 
(28) 

(29) 



f-l(x m ax,t + 1) = f+l(x max ,t) + SUC, 



(30) 



which ensures & constant material influx ji n — u(j)i n . On 
the other end of the simulation we have a free advection 
outflow boundary condition 



/ + i(0,t + l) = /_i(0,t)-«4(0 ) t), 



(31) 



ensuring an outflow of j out — u(f>(0). Both boundary con- 
ditions are implemented as the densities are advected. 
The boundary conditions are open and therefore the to- 
tal integral of the order-parameter is not conserved. To 



5 



calculate the derivatives of the concentration required for 
the calculation of [i in Eq. ([3]) at the boundaries, we de- 
fine the concentrations outside the lattice to be 

4>( x max +1,*) = 4>in , (32) 

0(-l,i) = 0(0, t) • (33) 
This fully defines our lattice Boltzmann method. 



IV. DIMENSIONAL ANALYSIS 

Our main purpose is to predict the final morphology 
formed by enslaved phase-separation fronts. In Fig. [T] 
we show the evolution of the phase-separation front and 
we see that only regular alternating domains are formed. 
Thus, the final morphology is fully characterized by the 
length of the domains formed as a function of the param- 
eters: 

A = \(a M ,as,b M ,bs,mM,ms,KM,Ks,CM,cs,u,(j) m ) . 

(34) 

However, not all of these parameters are independent. 
To reduce the number of free parameters we define the 
following nondimcnsional scales: 



X = 



-a s 



2tt V 2k s 



= _t_ = tm s a s 



"8p 



<Peq 



b s 



-a s 



(35) 



(36) 



(37) 



Here the length and time are scaled by their spinodal 
values, and the concentration is scaled by the positive 
equilibrium concentration of the separated material. In 
contrast with homogeneous quench spinodal decompo- 
sition and free-front propagation, the initial nondimcn- 
sional concentration (<j)i n ) does not play a meaningful roll 
in the scaling of the dynamics of phase-separation for en- 
slaved front systems. We therefore set (j>i n — for the 
definition of spinodal wavelength in Eq. (J7J) and spinodal 
time in Eq. ([8]) used for nondimensionalization. 

When we apply these nondimensionalizations to the 
equation of motion (jlip we obtain: 



a(x) 



1 / mix) 
- <l> 1 = 2^ Vx ' V * 



m s 



(38) 



as 



bs 



1 K(x) 2 , c(x) 



8ir 2 ks 



-as<t>eq 



This leaves us with the following dimensionless parame- 
ters: 



U = 

M = 
K 



2k s 



m M 
m s 

KM 



nms V (-a s ) 3 



, A 

C = 



as 
cm - c s 



D 



b_M 

bs 



a s 4>. 



(39) 



(40) 



The main effect of &a/ is to renormalize ajvf and cm, 
and the additional nonlinear effects can typically be ne- 
glected unless the concentration varies substantially in 
the mixing region. We therefore consider B an irrelevant 
nondimcnsional parameter and choose it to be zero. We 
choose C in such a way as to avoid bulk chemical poten- 
tial differences between the mixed and separated regions 
fi\j lk (x — > oo) = ^ h g lk (x — > — oo), and we obtain 



C = A$ in . 



(41) 



Additionally, we do not consider changes in the interfacial 
free energy cost across the front, and therefore set K — 
1. This results in a nondimensional domain wavelength 
which is a function of four nondimensional parameters: 



L = L(A, M,U,$ in ) 



(42) 



<:</ 



To investigate the dependence of L on U we first turn 
to numerical simulations. 



V. SIMULATION RESULTS 

Let us first consider a generic simulation to discuss the 
main features of morphology formation. The result of 
one simulation using this one-dimensional LBM imple- 
mentation is shown in Fig. [T] as a space-time plot of con- 
centration and chemical potential. It is clear from this 
plot that the dynamics of phase-separation is limited to 
the region around the front, and the chemical potential 
away from the front is flat. This simulation started with 
equal parts ^4-type and £?-type material in the mixed re- 
gion, and the A and B domains formed after the front 
passed are therefore of equal size. 

We shall now examine the domain formation using the 
middle .4-type domain in Fig. [T] as an example. Three 
key events in the formation of this domain are marked 
A, B, and C. At point A the domain is nucleated. The 
concentration plot shows a very narrow domain and the 
chemical potential plot shows the previous peak collaps- 
ing toward a flat profile. The domain there expands as 
it is pulled along with the front. To grow this domain, 
S-type material must be transported away. From the gra- 
dient on the chemical potential we see that most of this 
material is pulled across and deposited on the back of 
the domain, which advances the back interface of the do- 
main. Some of the rest is pushed ahead of the front, and 
a small amount builds up on the forming domain just be- 
hind the front. These last two depositions have the effect 
of making the forming domain increasingly unstable as 
the concentration at the front increasingly deviates from 
the equilibrium concentration. To maintain the currents 
of material away from the interface, the chemical poten- 
tial gradients must be maintained. Therefore the trough 
in the chemical potential deepens. 

At point B the ^4-type domain has just detached from 
the front as a new 23-type domain has been nucleated. 
The chemical potential very quickly becomes flat across 
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Lattice Site 

(a) Concentration (f>(x,t) 




Lattice Site 

(b)Chemical Potential fi(x,t) 



FIG. 1: (color online) Timelapse plots for an enslaved phase- 
separation front moving from right to left at nondimensional 
speed U = 0.001 leaving regular alternating domains in its 
wake. The first recorded profile is at the bottom of the graph 
with subsequent profiles shifted in the positive t/-axis. Each 
profile is separated by 200,000 iterations, or a nondimensional 
time difference AT = 12,500. The x-axis scale is in lattice 
sites. The upper graph shows the concentration profile over 
nearly two cycles of domain formation. The concentration 
ahead of the front can be seen to increase over time as the 
width of the two domains behind the front increases, until 
a new domain is nucleated. The lower graph is the total 
chemical potential for this system, showing that the gradient 
of the chemical potential is flat except for just ahead of the 
front and across the first domain behind the front. Additional 
simulation parameters are A — 1, M = 0.1, k = 2, cj>i n = 0, 
and s = 0.1. 



the newly detached A-type domain as it fully separates 
to the positive equilibrium concentration. However, the 
formation of a 23-type domain at the front causes A-type 
material to be deposited on the leading edge of the de- 
tached domain, thereby expanding it. The interface be- 
tween the A-type domain and the new S-type gets pulled 



along by the front, but is moving with about half the 
front's speed. The rate of deposition of material on this 
moving interface is nearly constant as revealed by the 
almost constant linear slope of the chemical potential. 

Point C marks the end of the formation of the domain 
by the enslaved phase-separation front as yet another do- 
main is nucleated. The chemical potential curve has be- 
come flat immediately adjacent to the domain, and all 
material currents to the interfaces have stopped. The do- 
main is now completely stationary. Because the domains 
are highly ordered and much larger than the interface 
width, they do not evolve further at any appreciable rate 

a 

We now qualitatively understand how the morphology 
is formed by a one-dimensional enslaved front. To make 
quantitative statements we need to measure the relevant 
quantities from the LBM simulations of our model. 



A. Measurement Methods 

The sizes of domains formed are found numerically by 
determining the length between zero-crossings of the con- 
centration. To find the zero-crossing point to sub-lattice 
precision we linearly interpolate between the two lattice 
sites on either side of the crossing. The first two domains 
after the control parameter front are still forming and ex- 
panding, therefore domain length measurements begin at 
the third zero-crossing behind the front and is an average 
of two or more domains. To ensure that there are enough 
fully- formed domains, we dynamically grow the simula- 
tion to always keep the number of interfaces behind the 
front greater than those used to calculate a meaningful 
average. The measurement is recorded once the uncer- 
tainty in the domain length is both less than one lattice 
spacing, and less than one thousandth the domain length 
itself. To more directly compare our enslaved front sys- 
tem to previous work, we double the domain size to find 
the domain wavelength. For asymmetric {4n n ^ 0) the 
domain size is multivalued and may lead to some con- 
fusion. In this paper, however, we only present domain 
wavelength measurements for symmetricly mixed mate- 
rial. 

In addition to the sizes of domains formed, we need to 
know the time-dependent concentration of the partially 
phase-separated material near the front. We cannot sim- 
ply measure the concentration at the front, as it would 
always be near zero due to the interface there. How- 
ever, since the chemical potential is continuous across 
the front we can invert the bulk chemical potential (de- 
fined as Eq. for k = 0) to find the relevant concen- 
tration. As shown in the timelapse plots of Fig. [TJ the 
chemical potential has a local extremum at the location 
of the control parameter front. The chemical potential 
value we use to find the near-front concentration is the 
interlattice extremum value of a polynomial fitted to the 
three most-extreme values near the front. 
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(a)Domain sizes as a function of u and k. 



10" 



II 



10 




10" 



U=u/u 



10" 



sp 



(b)Nondimonsionalized domain wavelengths. 

FIG. 2: (color online) Simulation results for domain wave- 
length A as a function of phase-separation front speed u 
for various interfacial energy costs. Results are shown for 
raw data (top) and after rescaling to nondimensional lengths 
and speeds (bottom). A single constant value of mobility 
rriM = m$ = 1/2 was used for these simulations. 



B. Domain Size as a Function of Front Speed 

Since it is the speed of the control parameter front 
that determines whether the phase-separation front is 
enslaved, and for many physical systems the front speed 
may vary over time, it seems natural to first observe how 
the front speed determines the size of domains. In addi- 
tion, the interfacial free energy cost k plays an important 
roll in determining homogeneous and free-front length 
scales. Therefore its effect on domain size will also be 
investigated. 



As shown in Fig. 2(a) domains get larger the slower 
the front becomes. Also evident is that the larger the 
domains become the smaller the effect from the interfa- 
cial free energy, and that for slow enough front speeds 
the effect would become negligible. However, we know 



from our nondimensionalization of the equation of mo- 
tion that all dependence on the interfacial energy cost 
should scale away completely, but only if the simulation 
is capable of faithfully reproducing the equation of mo- 
tion. In Fig. 2(b) the nondimensional rescaling of the top 
plot is shown, revealing an impressive data collapse onto 
a single curve. 

The data collapse provided by the nondimensionaliza- 
tion, though not unexpected, is very encouraging. We 
can be confident that the appropriate parameters for con- 
tinued investigation of this model are those revealed in 
Eq. (I40p by the nondimensionalization procedure. In the 
remainder of this paper, we investigate the dependence 
of the domain wavelength on the nondimensional param- 
eters given in Eq. (|42|) . We now analytically investigate a 
situation where the nondimensional front speed U is the 
only relevant nondimensional parameter. 



VI. ANALYTICAL SOLUTION FOR M=0 

To develop an analytical theory of domain formation 
we focus on the special case of M = 0, where there is no 
dynamics in the mixed material ahead of the front. This 
reduces the relevant parameters to just the enslaved front 
speed U. All other dimensionless parameters are irrele- 
vant since there is no dynamics ahead of the front. Due to 
its simplicity, this system lends itself especially well to an 
analytical approach. We will predict the resultant mor- 
phology (domain wavelengths) by explicitly solving the 
dynamics of phase-separation at the front. We accom- 
plish this in two parts: determine the time-dependence 
of the concentration just behind the front, and find the 
concentration at which a new domain will nucleate. By 
combining these two results we find the switching time 
t sw , and thus the domain wavelength: 



A = 2ut* 



(43) 



To numerically verify our analytical results we cannot 
simply set the mobility of the mixed region to zero. This 
would require that tm = 1/2, which would lead to an 
unstable simulation 15]. Instead, the M = case is im- 
plemented by placing the control parameter front at the 
positive boundary, and denning the off-lattice concentra- 
tion as (bin . 



Concentration at the Front 



To determine the time evolution of the concentration at 
the phase-separation front we examine the current of ma- 
terial into and out of the domain interface formed at the 
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FIG. 3: Real prin ciple branch of the Lambert W function for 
1/e < x < e [ly] . The Lambert W function appears in the 
solution for an enslaved front with M = &i„ = for U <^ 1. 




FIG. 4: (color online) Comparison of the simulation results 
for the time evolution of the concentration at the phase- 
separation front for fast and slow front speed with <j>i n — 0. 
Note that the fast phase-separation front speed simulation 
shows marked deviation from the theoretical curve, as ex- 
pected. Also note that the concentration at which the domain 
switches (cj> aw ) does not seem to depend on front speed. Inset 
shows the full theoretical curve, though concentration should 
never go below <j> sp as that is unconditionally unstable and 
would immediately undergo spinodal decomposition. While 
not essential to this analysis, zero concentration occurs for 
TU 2 = (ln(16) - 2)/tt 2 » 0.078. 



front. We first assume the interface to be sharp by com- 
parison to the size of the domain forming at the phase- 
separation front. From the definition of the nondimen- 
sional speed U in Eq. (|59")) we can see that U oc s/k. 
Thus a vanishing cost of interfacial free energy (k — ► 0) 
is equivalent to assuming the front speed is much slower 
than the spinodal velocity, i.e. [/ « 1. This assump- 
tion of a slow front speed is not especially limiting to 
our analysis as we are primarily interested in the for- 



mation of large stable structures; unlike those found for 
front speeds approaching the free-front speed u* which 
are small and rapidly coarsen ■ 

The buildup rate of material in a region is 



J 



dt4> dx 



V-jcfc , 



(44) 



where the current j at a point is determined by the gra- 
dient of the chemical potential: 



j(t) = -mVM(t) 



(45) 



We're interested in the region of the interface forming 
directly behind the front at position x = xq. Since that 
region is vanishingly small, the rate of material deposited 
there can be written as a difference of the chemical po- 
tential gradients on either side of the front 



J 



(mAfVfiM - m s V^ s )S{x - x ) dx 
m M V//M - m s Vns ■ (46) 



Since the mobility ahead of the front is zero, the time 
dependent rate of deposition for the interface forming at 
the front is: 



J(t) = -msVAt(t) • 



(47) 



The gradient of the chemical potential across the entire 
forming domain is nearly constant, as can be seen in 
Fig. [TJ This means that material pulled behind the front 
is transported all the way across the forming domain and 
deposited on the second domain. Any deviation from a 
linear profile of the chemical potential is of order dt4>, as 
can be seen from the Cahn-Hilliard equation ||3J| and the 
fact that (f>(x, t) is nearly stationary inside the forming 
domain. 

The value of the chemical potential at the phase- 
separation front /io changes with time, but the chemical 
potential at the other end of the domain is zero. We can 
therefore write the chemical potential gradient as 



i{t) 



(48) 



where t is the width of the domain forming at the phase- 
separation front. 

The interface at the phase-separation front moves with 
the front, expanding the domain forming behind it. This 
requires that the rate of material deposited into the 
phase-separation front interface be equal to the differ- 
ence in concentration across the interface times the front 
speed 



J(t)=u [<j)o{t) - (j>i 



(49) 



where u is the (dimensional) front speed, 4>o{t) is the 
concentration at the phase-separation front, and is 
the mixed-material concentration. 
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The M = condition means that none of the 
wrong-type material can be pushed ahead of the phase- 
separation front; instead it must be transfered across and 
deposited behind the forming domain. This results in the 
expansion of the second domain by moving the interface 
at a speed ui < u which trails the phase-separation front. 
Taking this into account, and using an argument similar 
to the one used to derive Eq. (|49p. we can write another 
expression for the rate of material deposition behind the 
front in terms of material removed and transported to 
the second interface. We obtain 



J(t) = (u-ui(t))2<t> eq 



(50) 



where 2fa q is the change in concentration across the fully 
formed interface, and <f> eq is the equilibrium concentra- 
tion given by Eq. (j4]). For now we'll continue to write 
fa q , and eliminate the variable later through nondimen- 
sionalization. 

The width of the first domain is found by integrating 
over time the relative speeds of the interfaces at both 
ends of the domain 



rt ft 
£(t) = I ui(t)dt = u I 11 



<Po(t) - fa. 



20, 



eq 



dt , (51) 



where we eliminated uj by combining Eq. (|49|) and 
Eq. dSO]). We then combine Eq. (gTJ) with Eq. gl]), use 
the definition of Vfj,(t) from Eq. (|4"8j). then subtitute the 
definition of the chemical potential for our model, to find 
the following alternative expression for the forming do- 
main's width: 



i{t) 



u fa(t) - fa 



m s asfajt) + bg^pft) 
u fa(t) - fa. 



(52) 

Equating Eq. (fS"Tj) with Eq. (|52p we find the integral equa- 
tion 



1 



fa (t) - fa 



dt 



ms asfajt) + b s fa](t) 
u 2 (t> (t) - fa, 



for the concentration at the phase-separation front as a 
function of time. Differentiation of Eq. ([55]) with respect 
to time yields 



1+- 



2fa q 

Therefore 
d t fa = 



„2 



9Q 




-u 2 (4> - m ) 2 (2cj) eq + fa - <j> in ) 
2ms4>e q (bs^o (20o - 30, n ) - as fan) 



(55) 



where the time dependence of fa is implied. We nondi- 
mensionalize this using Eqs. (|35ti37[) . and obtain 



d T <S>o =k z U 



2tt2 ($ - $ m ) 2 (2 + $ - * m )(3<i> 2 „ - 1) 



2$g - 3*g$ in - 



(56) 
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FIG. 5: (color online) Free energy tangent construction for the 
switching condition for an enslaved phase-separation front. 
The nucleation kernel has concentration <f>i n and domain type 
switching will occur when the concentration just behind the 
front reaches <f> sw . Also marked are the equilibrium (4> eq ) and 



spinodal 



concentrations. 



which, for general U and $i„, can be solved numerically. 

For the = case of symmetric inflow concentra- 
tion, Eq. (p)6"|) simplifies to 



<9 T $ = -ir 2 U 2 



1 1 

+ 2 



(57) 



The initial concentration must be at equilibrium 
($o(0) = 1), allowing us to find the analytical solution 



<!V'f; : 2 + 21F j -I eX pQ^ 2 C/ 2 T-i 



(58) 



where W is the principle branch of the Lambert W func- 
tion. The Lambert W function is the solution to the 
equation x = Wix) exp(W(x)) for complex x and has in- 
finitely many branches. The principle real branch used 
for this solution is plotted in Fig. [3J and more information 
on the Lambert W function can be found in (l6j . 

As explained earlier, our simulations are able to track 
and record the concentration at the front by inverting the 
chemical potential. In Fig. 2] we show that the analytic 
solution is in excellent agreement with our simulation 
results for slow phase-separation front speeds. For faster 
fronts there is a deviation which is consistent with our 
assumption of small U. 



B. Switching Condition 



10 



— Switching Theory 
O Simulation 




'sp 

-(b 



FIG. 6: (color online) Verification of the switching condi- 
tion for various mixed material concentrations. The solid line 
shows the analytic solution for the value of (j> behind the front 
at which switching occurs. The circles are switching concen- 
trations recorded from LBM simulations. The simulation pa- 
rameters were M = 0, and U = 6.366 x 10 -3 . 



We will now need to determine when a new domain will 
be formed. At first glance one might expect that the sys- 
tem will switch when the concentration reaches its spin- 
odal value. As we saw in Fig. [4] however, the switch 
typically occurs much earlier. Consider an ,4-rich do- 
main forming behind the front. Overtaken material can 
either extend the existing domain by transporting away 
excess £?-type material, or nucleate a new £>-rich domain; 
the system will choose the option that minimizes the to- 
tal free energy. We now observe that there is already an 
interface between the domain forming behind the front 
and the material the front is overtaking. Because the 
interface required to nucleate a new domain is already 
in place, we only need to consider the bulk free energy 
for the switching mechanism. The inflowing material will 
nucleate a new domain, even without fluctuations, when 
such a domain formation continuously minimizes the bulk 
free energy. 

This argument is sketched in Fig. [5] The concentra- 
tion of the potential nucleus of B-type material forms one 
point on the free energy curve. The other point on the 
free energy curve is the concentration of phase-separated 
A material just behind the front. If the line-segment in- 
tersects the free energy curve at some other point, the free 
energy can be lowered by dissolving the potential type- 
B nucleus, thus diluting the domain of type-„4 material. 
Nucleation will occur if the line-segment connecting these 
two points on the free energy curve is below the free en- 
ergy curve. Thus the system will always switch before 
the spinodal concentration is reached. 

In our case the mixed material ahead of the control- 
parameter front has a vanishingly small mobility (M = 
0) , so the effective concentration of material ahead of the 
front is simply the concentration of the mixed material. 



With that point firmly fixed on the free energy curve, we 
can predict the critical concentration at which domain 
switching will occur. 

We denote the switching concentration — the concen- 
tration at which domain switching due to nucleation 
occurs — as <p sw . We just need to find the line which lies 
tangent to <j) sw on the free energy curve and crosses the 
nucleation kernel concentration <^ n . The slope of this 
line is the slope of the free energy curve evaluated at the 
switching concentration. This is also the definition of the 
chemical potential at the switching concentration. 



F( 



,)-F{4> in ) dF 



(59) 



This construction is illustrated in Fig. and results in 
the closed-form solution for the switching concentration: 



± 



(60) 



We verified Eq. (jST))) using our LBM simulations by 
tracking the peak of the chemical potential near the front 
during a complete half-cycle and recording the maxi- 
mum value of the peak for different inflow concentra- 
tions <pi n . This is then inverted to find the concentra- 
tion as a function of the chemical potential. The con- 
centration at the maximum value of the chemical poten- 
tial peak is the switching concentration. The results of 
this test are shown in Fig. [6l and demonstrate excellent 
agreement with the theoretical prediction in the range 
of — (j) sp < (f>i n < 4> sp . Outside this range the system is 
meta-stable and we do not find alternating domains. 



Domain Size 



The expression for the switching concentration in 
Eq. (|60|) is trivially nondimensionalized 




as is the expected domain wavelength 
L = 2UT SW . 



(61) 



(62) 



Here T sw is the nondimensional switching time defined as 
the time when the concentration at the front reaches the 
switching concentration. Plugging the switching concen- 
tration of Eq. (|61[) for symmetric mixed material inflow 
(&in = 0) into the solution of the differential equation 
for the concentration at the front, presented in Eq. (|58|) . 
yields 



2 _ 2 



2W 



(63) 
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FIG. 7: (color online) Nondimensional domain wavelengths 
for M = from lattice Boltzmann method simulations 
agree very well with the analytical prediction for low phase- 
separation nondimensional front speed U. Note that the the- 
ory contains no adjustable parameters. 



which can be inverted to solve for the switching time: 

2 (%/6 + 61n (2- y/2/3) - s) 
Ts W = ^ ^ ■ (64) 

With this value for the switching time we can, at long 
last, explicitly solve for the domain wavelength 



L(U) = 



4(V6 + 61n(2 



(65) 



3tt 2 [/ 

as a function of phase-separation front speed for sym- 
metric inflow {4>in — 0). This solution contains no free 
parameters and only relies on the front speed being much 
smaller than the spinodal speed. 

We verify Eq. (|65|) by comparing it to data taken 
from the LBM simulations. The results are shown in 
Fig. [7] We see that the simulation results agree ex- 
cellently with our theoretical parameter-free prediction 
for nondimensional phase-separation front speeds slower 
than U ~ 1CT 3 . 

The domain wavelengths for asymmetric inflow (<&i ra ^ 
0) can similarly be obtained by numerically solving 
Eq. (|56|) for the switching time as a function of <fii n . 



VII. SOME PROPERTIES OF MORE GENERAL 
CASES 



We have successfully developed an analytical expres- 
sion for the morphology formed by an enslaved phase- 
separation front in the special case where there is no 
mixed material dynamics and the mixed material is sym- 
metrically composed. We accomplished this by con- 
structing a differential equation from the analysis of cur- 
rents which drive the dynamics of phase-separation at 
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FIG. 8: (color online) Dimensionless domain wavelength L 
as a function of mobility ratio M as measured from LBM 
simulations. The inset plot shows the same data over larger 
mobility ratio range. This shows how the family of curves 
found in the L vs. U plane depends on M. For the red 
(online) circles, the phase-separation front speed was kept a 
constant at U = 6.366 x 10~ 3 , while k, ms and tum were 
varied. The blue (online) square is from a simulation of the 
special case M = at the same nondimensional speed. Note 
this curve is not monotonic. 



the front. We could only solve this differential equation 
analytically after the simplification allowed by assum- 
ing symmetric mixed material, however the more general 
case can be tackled numerically. Along these same lines, 
one should be able to construct an even more general dif- 
ferential equation which includes dynamics ahead of the 
control-parameter front. However, even without such an 
in-depth analysis we can, at least qualitatively, under- 
stand the effects of mixed material mobility on domain 
formation. 

We "turn-on" the dynamics ahead of the control pa- 
rameter front by setting the nondimensional mobility M 
to some positive value. This results in a family of curves 
in the L vs. U plane, each of which are similar in shape to 
Fig. [2 To see how this affects the size of domains formed 
we can choose some slow front speed and perform a se- 
ries of simulations with varying M . The results of such 
a series of simulations is shown in Fig. [5J We observe a 
rapid but continuous reduction in domain wavelength as 
the nondimensional mobility is increased from zero, until 
it reaches a minimum. We should mention that data of 
two different kinds of simulation are shown in the graph. 
The simulation for M — where the front is at the inflow 
boundary, and M > where the front is firmly inside the 
lattice. There is excellent agreement between the meth- 
ods. 

After the sharp reduction levels off, L increases with 
M and eventually results in the formation of very large 
domains. This nonmonotonic behavior is discussed be- 
low, followed by an explanation of the effects of the final 
free nondimensional parameter A. 
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A. Small Mixed Mobility (M<1) 

If the dynamics of material ahead of the front are very 
slow by comparison to the separated region, the currents 
are almost identical to the case of no mixed-material dy- 
namics. Such a case could be physically realized for a sys- 
tem where the mobility strongly increases with tempera- 
ture and we have a lower critical point. In this case, very 
little material is allowed to build up ahead of the front. 
Currents ahead of the front are then vanishingly small. 
What does change, due to this small buildup ahead of the 
front, is that the front is now overtaking a slightly larger 
nucleus of material. This slightly larger nucleus will make 
domain type switching more energetically favorable. 

We can see the effect a small increase in nucleus vol- 
ume has on the switching condition by re-examination of 
Fig. [5j A slightly larger nucleation kernel has the effect 
of shifting <pi n towards the opposite material type. This 
results in a shallower tangent line construction which con- 
tacts the free energy curve closer to the equilibrium con- 
centration. This means that the concentration of ma- 
terial just behind the front cannot be diluted as much 
before domain type switching occurs, and earlier domain 
type switching results in smaller domains. 

B. Large Mixed Mobility (M>1) 



This is a case that is easily realizable in real systems. 
For strongly phase-separated regions the mobility can be- 
come very small compared to the mixed material. As we 
saw from the results of the simulations shown in Fig. [5J 
for large M the domain wavelengths become very large. 
This happens because diffusion is rapid in the mixed ma- 
terial region, and large amounts of material can build 
up ahead of the front. When the domain at the front 
switches type, this buildup floods into the new domain, 
thereby making it very large. In the limit of large M, 
the domain is predominantly formed from this material. 
Since this is effectively the same as halting the dynamics 
behind the front, we can understand this further by ex- 
amining just the mixed material dynamics ahead of the 
front. 

In the reference frame of the front, material pushed 
ahead of the front is governed by the drift-diffusion equa- 
tion 

d t (f> + V(u4>) = VtoV0 , (66) 

which is defined for x > xq. Here we have neglected 
the interfacial term and absorb the a parameter into the 
mobility. The boundary condition at the front at x = xq 
is 

J (t) = u[4>(x , t) - 4> eq ] , (67) 

which is the current of material rejected by the forming 
domain, analogous to the result of Eq. [)49p. This bound- 
ary condition results in a buildup of material which will 




M 

FIG. 9: (color online) Domain length scaling is universal for 
changes in front speed and mobility ratios in the slow front 
speed limit. The upper graph shows that the nondimensional 
domain wavelengths L and front speed U both vary by an 
order of magnitude, yet retain similar curve shapes. 

eventually lead to switching. These two equations are 
sufficient to define the material dynamics ahead of the 
front. We now nondimensionalize these equations to re- 
cover their parameter less versions. 

We introduce the general length and time scales 

x = Xx , and t = Tt . (68) 

It should be noted that X and T are not the same as 
the spinodal wavelength and time we used in the earlier 
nondimcnsionalization. The drift diffusion equation then 
becomes 

^+^v(^)=lv(r,iw) , (69) 
which we rewrite as 

a^ + v(<ty) = v(^wj . (70) 

For the boundary condition 

J(t) = -u(4>(x ,t) - (f> eq ) , (71) 
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we obtain the nondimensionalized boundary current: 

J(t) = <u(<f>(x ,t) ~<j> eq ) . (72) 
We now choose the length and time scales such that 



l = u= —u , and 1 = 



(73) 



This results in the parameterless equation of motion 

d t <j) + = V 2 , (74) 
and boundary condition 



J = (4>(x ,t) - (f> eq ) 



(75) 



While we do not have the analytical solution of this 
differential equation, we know that the solution exists 
and will result in a nondimensional switching time t sw . 
We use this to replace the dimensional switching time 
appearing in Eq. (|4*3")) , which gives the domain wavelength 



A = 2uT t s 



(76) 



We now evaluate T in terms of m and u from Eq. ([73 



T 



m 
1 1 ^ 



This uncovers the expected domain wavelength 



A 



in 



m - 
2 t sw oc 
u u 



(77) 



(78) 



We can re-express this result in terms of the earlier nondi- 
mensionalization of Eqs. (|35H37p as 



M 
Loc-. 



(79) 



Thus we uncover the expected scaling behavior of the 
domain wavelength for the simultaneous limit of small U 
and large M. 

Interestingly, the scaling behavior of 1-D coarsening of 
a homogeneous quench is logarithmic in time, therefore 
enslaved phase-separation fronts can build large domains 
much more effectively than the phase-ordering which fol- 
lows spinodal decomposition. In principle, the speed of 
the front can be controlled, and thus the production of 
highly ordered, depth-dependant structures becomes pos- 
sible. 

We have already seen from the analytical result for 
M = in the limit as U — * that the domain wavelength 
is proportional to the inverse of the front speed (L oc 
U^ 1 ) which matches the previous scaling argument for 
the opposite case of M — * oo. The scaling argument 
reveals that LU vs M will result in universal scaling in 
the limit of very slow fronts. We test this by performing 
simulations on a selected range of sufficiently slow front 
speeds U < 10 -3 , and a range of mobilities. The results 
of these simulations are shown in Fig. [§] revealing the 
large difference in domain wavelengths in the U vs M 
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FIG. 10: (color online) Switching concentration for large 
nondimensional mobility M as a function of control parame- 
ter height A. The predicted curves are calculated by assuming 
perfect continuity of the chemical potential across the front, 
and using this to predict the size of the nucleation kernel in 
the switching condition. The dashed line includes the addi- 
tion of a current j which is caused by the dynamics ahead 
of the front, and is zero in the M = case. The simulation 
used a nondimensional mobility ratio M = 5, and front speed 

c/ = io- 4 . 



plot which then collapse reasonably well in the LU vs M 
plot. 

To improve the collapse, even slower simulations would 
be required, however this becomes too computationally 
expensive. For the measurements presented here we re- 
quired simulations of more than 10 5 lattice sites simu- 
lated for more than 10 8 iterations. Such simulations take 
about one week on a Xeon powered Linux workstation. 



C. A and M at the Front Boundary 



The relative quench depth A is the last relevant parame- 
ter which needs to be considered in determining morphol- 
ogy formed by slow enslaved phase-separation fronts in 
ID. The control parameter in either a purely mixed (Am) 
or purely separating (As) system is usually not a func- 
tion of position and can be absorbed into the mobility. 
In fact, if we were to consider the mixing region ahead of 
the front as being entirely detached from the separating 
region behind the front, we would be able to completely 
nondimensionalize away all free parameters, resulting in 
separate scale-invariant equations of motion. However, 
due to the shared boundary at the front linking the two 
regions, the change of the control parameter A across the 
front cannot be ignored. This means that changing the 
control parameter ratio A has different consequences for 
the resulting morphology than does changing the mobil- 
ity ratio M. 



14 



The first consequence of changing the order parameter 
A has for the dynamics of phase-separation is for the ma- 
terial which gets pushed ahead of the control parameter 
front. Recall that the chemical potential is continuous 
across the front (/ig = /im)- Ignoring for the moment 
the interfacial contribution to the chemical potential, we 
can write the continuity of the the chemical potential as 

0s - 4>s = Mm , 

=► <t>M = - (4% - 4>s) ■ (80) 

This means that the concentration just ahead of the front 
and the concentration just behind the front are tied to- 
gether. This is important because the switching con- 
dition, even for very large M, is partly dependent on 
the concentration of material just ahead of the front. 
For A <C 1, the concentration ahead of the front be- 
comes larger, causing the switching concentration to ap- 
proach the equilibrium concentration, resulting in very 
early switching. On the other hand, for A ^> 1 the 
concentration ahead of the front becomes close to the 
symmetrically mixed concentration, resulting in switch- 
ing taking longer to occur. 

If we simply use 4>m in Eq. (|80p as an effective <^ n in 
the earlier derived switching condition of Eq. (jHO)) , we 
discover a new switching condition for large M which 
is dependent only on A. This predicted switching con- 
centration is plotted as the solid line in Fig. along 
with switching concentrations measured from our LBM 
simulations. Since the switching condition assumed no 
dynamics ahead of the front it is not surprising that this 
prediction is inadequate for the opposite scenerio. If we 
take into account the additional current j across the front 
due to the dynamics ahead of the front, and account for 
this in the effective <pi„ in the switching condition 

<t>t f = j(<f>%-<t>s-j) , (81) 

we can measure this current at one value of A and use 
it to fit the rest of the curve. This new prediction is 
the dashed line in Fig. [101 with the A = 1 switching 
concentration used as the fitting value. This new predic- 
tion agrees well, although not perfectly, with the simu- 
lated values, showing that the current across the front at 
switching seems to scale as the inverse of A. 

VIII. OUTLOOK 

Phase-separation front driven morphology is a rich and 
complex subject which has significant potential for new 
research and applications. In this paper we focused on 
the simplest case: a one-dimensional binary mixture. For 
abrupt control-parameters we found that a regular al- 
ternating morphology is formed. We have shown that, 
complex though the subject may be, certain important 
aspects of enslaved phase-separation fronts can be under- 
stood well enough to theoretically predict the morphol- 
ogy 



Using our model of a binary mixture with an abrupt 
control parameter front, we were able to determine the ef- 
fect of all relevant free parameters on the size of domains 
formed in a ID system. We verified our predictions by 
comparing them to LBM simulations of our model. 

We were able to reduce the number of relevant param- 
eters from twelve to four. We achieved this by rescaling 
length and time with the spinodal wavelength and time 
and the concentration by the equilibrium value of the 
order parameter, and by making some simple physical 
assumptions. 

Analytical solutions are few and far between for the 
theory of morphology formation. Choosing the nondi- 
mensional front speed as the primary parameter of inter- 
est, we consider the radical simplification of a front mov- 
ing into material which has no mixed-material dynam- 
ics. This fixes the other free parameters and allows us to 
analytically determine the dynamics of phase-separation 
induced by a slow front. 

We first find the time-dependence of the concentra- 
tion at the front. We then determine the mechanism for 
domain switching. While typical nucleation arguments 
are heavily dependent on the interfacial free energy our 
switching condition was found to depend purely on the 
bulk free energy. 

This gives us an analytical prediction of the formed 
which is in remarkable agreement with the numerical re- 
sults. Additionally, this solution functioned as a stepping 
stone for understanding the effects of the remaining pa- 
rameters. 

By allowing dynamics in the material ahead of the 
front we discovered a nonmonotonic dependence of the 
domain wavelength on the nondimcnsional mobility. Al- 
lowing a small amount of material to be pushed ahead of 
the front actually decreased domain wavelength, whereas 
allowing a large buildup of material ahead of the front 
resulted in huge increases in domain wavelength. The 
former is understood by an extension of the switching 
condition, and the later is explained by a simple scal- 
ing argument; when the mobility in the mixing region is 
much greater than separating region, the domain wave- 
length formed by the front should scale with the same de- 
pendence as the length scale in the mixing region. This 
argument resulted in the discovery that, for very slow 
front speeds, there is a single universal scaling curve of 
the nondimensional domain wavelength as a function of 
the nondimcnsional mobility and front speed. 

Domain type switching is very well understood for the 
case of no mixed material dynamics, and we've explained 
the effect of very small mixed material mobility. The 
effect of large mixed material mobility on the switching 
condition is uncovered in the analysis of the quench depth 
A of the enslaved front. Modification of the switching 
condition is accomplished by including the continuity of 
the chemical potential, and a current j which is inversely 
proportional to A. 

The effect on one-dimensional morphology formation 
of three of the relevant nondimensional parameters is now 
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well understood. We now briefly mention several research 
avenues that lead on from the results presented in this 
paper. 

We are currently examining higher-dimensional en- 
slaved control parameter fronts. Two-dimensional fronts 
open the possibility for the formation of stripe morpholo- 
gies which are perpendicular or at an angle to the control 
parameter front. Additionally we find other exotic mor- 
phologies like polka-dot lattices and ovoid domains. Also, 
dimensionalities of two and higher allow hydrodynamics 
to play a roll in the phase-separation dynamics. In three- 
dimensions, enslaved fronts can induce a rich family of 
morphologies depending on the properties of the material 
and the front. Perhaps the most intriguing possibility is 
to be able to control the material and front well enough 
to reliably switch morphologies during front traversal, 
resulting in highly-ordered, highly-inhomogeneous com- 



posite materials. 

What also remains is to examine more general en- 
slaved fronts. The most straightforward extension is to 
attempt to develop an analytical solution to the abrupt 
one-dimensional enslaved front for all three free parame- 
ters U, A, and M. A careful accounting of the additional 
currents cased by dynamics ahead of the front should al- 
low an extension of the differential equation describing 
the dynamics at the front, presented in Eq. (f56|) . for this 
more general abrupt one-dimensional enslaved front. 

Then we will relax the condition of an abrupt control 
parameter front to obtain an extended front. We ex- 
pect that there will be little change while the width of 
an extended control parameter front is on the order of 
the interface width, but extended fronts which approach 
the size of the forming domain may result in new and 
interesting phase-separation dynamics, even in ID. 
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